Creep and fluidity of a real granular packing near jamming 
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We study the internal dynamical processes taking place in a granular packing below yield stress. At 
all packing fractions and down to vanishingly low applied shear, a logarithmic creep is observed. The 
experiments are analyzed using a visco-elastic model which introduces an internal, time dependent, 
"fluidity" variable. For all experiments, the creep dynamics can be rescaled onto a unique curve 
which displays jamming at the random-close-packing limit. At each packing fraction, we measure 
a stress corresponding to the onset of internal granular reorganisation and a slowing down of the 
creep dynamics before the final yield. 

PACS numbers: 47.57.Gc, 83.80.Fg, 65.60.+a 



For granular matter, it is currently accepted that a 
quasi-static limit exists as for grains of macroscopic size, 
thermally activated processes can be ignored. At low 
shear rate, mechanical properties of granular packing are 
usually described by rate independant constitutive rela- 
tions [l|. However, there are compelling experimental 
evidences that this limit is just a short- time approxima- 
tion and time dependent processes are significant on the 
long run [H. Many numerical simulations based on soft 
interparticular interactions [sl have brought to the front 
the idea of an "universal" jamming transition scenario 
based on a mechanical rigidity threshold separating solid 
and fluid behavior (see a recent review and refs in ([sl)). 
However for real grains, interparticle solid friction was 
shown to affect the rigidity onset and stabilize packing 
at compacity below the random close packing limit [iq| . 
In this case, experiments have pointed out the central im- 
portance of nanometric scales where humidity [ll[, con- 
tact plasticity, tiny thermal variations [l2| or weak me- 
chanical external noise [l3|, do impact significantly the 
macroscopic dynamics and the rheology. Since the orig- 
inal theoretical propositions of " soft glassy rheology" or 
"'Shear Transformation Zones" (Hf various models have 
tried to capture the complex energy reorganization dy- 
namics taking place in amorphous solids or yield-stress 
fluids, in relation with their constitutive rheological laws. 
Microscopically, the emergence of plasticity is often ex- 
plained via a simple picture where localized elastic insta- 
bilities release irreversibly long range elastic constraints 
[sl which may be organized spatially as shear driven 
avalanches 0. Macroscopically, to account for the com- 
plex phenomenology, an internal time-dependent variable 
called "fluidity" is often introduced @ to describe the 
rate of stress relaxation. In this paper, we study the 
creeping dynamics of granular packing under constant 
shear stres, below the Coulomb limit. In order to reveal 
the internal relaxation processes and make connection 
with the behavior of a large class of yield stress fluids [t^ , 
we propose a quantitative analysis using a phenomeno- 
logical model based on the "fluidity" concept. 



2000 4000 
Times : (s) 



FIG. 1: (a) Schematics of the shear cell. (M): motor, (S): tor- 
sion spring, (T): torque probe, (D): induction distance probe, 
(A): transversal arm, (Q): flowmeter. Pi: differential pressure 
probes, (b) Creep deformation : Aj(t) under constant shear 
at packing fraction = 0.60. 

All the mechanical tests are performed at well defined 
packing fractions (j) (see figH]). To achieve this goal, the 
set-up is designed as an air fluidized bed. The container 
is a plastic cylinder of inner diameter D = 10cm closed 
at its bottom by a metal grid supported by a honeycomb 
grid. Pressurized air is introduced in an admission cham- 
ber below the grid at a controlled over-pressure AP. We 
use glass beads of density p = 2500kg /m^ and mean di- 
ameter d = 200 fim (R.M.S. polydispersity Ad = 30 fim). 
A mass M of grains is poured into the container such that 
the typical packing height is L = 10cm. Using a flow rate 
just above the fluidization value, we obtain after stop- 
page, an initially loose granular packing at a compacity 
(j) ^ 0.56. Then, by successive tapping on the container 
side, the packing fraction can be increased up to the de- 
sired packing fraction (maximal value (f) = 0.625). Note 
that in the present report, the relative humidity is kept 
at 35 ±5% and we work in a "quiet" environment charac- 
terized by a backgound noise quantified by placing in the 
packing an accelerometer (RMS acceleration 0.02ms~^ ). 
Before each mechanical measurement, the packing frac- 
tion value is evaluated from a linear fit between pressure 
drop AP and flow rate Q : AP = AQ, stemming from 
Darcy's law [13] • The relation between permeability K 
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FIG. 2: . (a) Shear modulus (^(0) under gravity confinement. Straight line y — Gox^ with Go = 2, 28.10^ Pa. (b) Response to 
a stress ramp : shear stress cr as a function time for /3 = 0.00104 rd/s and a packing fraction (j) = 0.6. First rupture stress ar 
and maximal stress ctm are displayed as horizontal dotted lines, (c) For the same rotation rate, dr, (Jm and the ratio as a 
function of packing fraction. 



and packing fraction was calibrated by preliminar series 
of experiments and a Carman-Kozeny relation was ob- 
tained : K{(l)) = Aii^d^, with A = 1/165. Conse- 
quently, for a mass M of grains poured in the cylinder, 
the packing fraction is obtained through the relation : 
(j) = V XT^D^p -) where 7? = 1, 85x10 — 6 Pa.s is the air 
viscosity. To shear the granular packing, we use a four- 
blade vane in stainless steel of height Hq = 2.54cm and 
diameter 2Ro = 2. Mem (see fig jl^), introduced at a depth 
h = 5cm below the surface prior to the initial fluidiza- 
tion process. This procedure creates reproducible initial 
conditions. Shear stress is applied through the vane (see 
fig CI connected axially to a torque probe (T) itself cou- 
pled to a brushless motor (M) via a torsion spring (S). 
The vane rotation angle a is monitored via a transversal 
arm (A) whose rotation is followed by a displacement in- 
duction probe (D). The motor rotation angle (3 is imposed 
with a 27r/ 10000 precision. Torque and displacement sig- 
nals as well as the motor command are connected to a 
Labview controler board. The last one is programmed 
to impose a motor rotation rate or a fixed stress using a 
feedback loop on the torque signal. In the following, we 
ignore the stress and strain spatial distribution due to 
the Couette cell geometry and define only average values 
obtained from mesurements of angular rotation a and 
torque T. The mean packing deformation 7 is defined 
as 7 = and the mean shear stress is a = 2'kR^Hq ' 

On figlTJb), we display three examples of creep curves 
A7(t) = 7(t) — 7(0) obtained at fixed compacity and 
shear stress values a. 

Elastic response - To obtain the elastic response of the 
packing initially prepared at a given packing fraction, 
stress cycles were performed corresponding to sinusoidal 
deformations of small amplitudes around 10 The cy- 
cles were carried out under constant mean confining pres- 
sure (hydrostatic loading). In the short time limit, the 
response is essentially reversible (see inset of fig.(j2^). 
The effective elastic shear modulus increases with pack- 
ing fraction almost linearly : Geff = G^cj) (see line on 



figEti), with a value of Gq = 2,28.10^Pa. Interestingly 
this simple result is consistent with a mean-field Hertz 
elasticity theory (seefl^ and refs inside) where under a 
confining pressure Pq = p(t>gh^ the shear modulus scales 

as Geff cx EQ{(j)Z)^^^ (^) ' ^^^^^ Z is a constant 
mean contact number and Eq the material Young's mod- 
ulus, meaning that, in the range of density explored, the 
contact density at the origin of the (j)Z term varies lin- 
early with (j) . 

Response to a stress ramp - To identify the the maxi- 
mal stress supported by the packing before yield, shear 
stress ramps were applied at a constant motor rota- 
tion rate (/3), using the softest spring constant available 
{k = 2,45.10"^A/'m/r(i). On figure [2fa), the stress re- 
sponse at (/) = 0.60 is displayed as a function of time. At 
first, a linear increase of the stress is observed with a slope 
corresponding to the spring constant. Then, at a given 
stress level <j^, we observe the emergence of well-marked 
and sudden granular material reorganizations (see top- 
view on figlSfb)) in the form of rather equidistant events 
corresponding to stress drops and large plastic deforma- 
tions (^7=10~^ — 10~^). We define this stress value as 
the "first rupture stress" : <j^. However, stress can still 
be increased but undergoes series of partial rupture, up 
to a maximal value <jm followed by a large stress jump. 
Then, a subsequent stick-slip dynamics is observed. In- 
terestingly, such a structured fluctuation regime was also 
observed by Albert et al.[l5| for the drag force on an in- 
truder near jamming. In our case, the maximal stress 
value corresponds to a Coulomb yield criterion as we 
verified that its value increases linearly with the confin- 
ing pressure. Such experiments were performed varying 
packing fraction and ramp velocities and we only dis- 
play here stresses obtained at the slowest driving veloc- 
ity where the values are quite insensitive to the rotation 
rate. On figlS^b), the rupture, maximal and dynamical 
stresses are displayed as a function of packing fraction for 
a rotation rate (/3 = O.OOlOArd/s). The values increase 
strongly with packing fraction. Interestingly these stress 
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values are strongly related since their ratio stay constant 
with packing fraction: X = aM/o'r = 3.5 ±0.2 (see upper 
part of fig.([2)D). 
Creep Bow - This part is concerned with the creep 
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FIG. 3: Creep experience, (a) plot of 7 ^ — 7q ^ as a function 
of time for creep experiments performed a constant packing 
fraction (f) = 0.6 for various shear stresses, the straight lines 
are linear fits : y = Ct . (b) Values of the fitted slopes C as 
a function of G{(j))/a for three values of packing fractions. 

response of the granular packing under constant shear 
stress. The packing fraction is varied between 0.56 and 
0.625. The procedure consists of two steps. First, a 
monotonic loading up to the desired stress value. This 
initial step is fast with respect to the creep dynamics, 
typically less than 200s. We verified that this loading 
time is much smaller than the inverse of the initial shear 
rate. Second, a phase of constant applied shear stress is 
obtained by a feed-back procedure where stress is main- 
tained at a constant value within a range less than 1 
%. The onset of feed-back defines the initial time t = 0. 
If stress leaves the assigned range, a command is sent 
to the motor to rotate the torsion spring and to adjust 
the torque accordingly. From time to time, as a con- 
sequence of this fast motor rotation, an acceleration of 
the strain rate is observed, followed by a decay down 
to the value before the jump. These dynamical phases 
lasting less than 20s were replaced by a linear interpo- 
lation of the shear rate so that they do not artificially 
perturb the subsequent analysis. For all the experiments, 
we observe a slow increase of the deformation Aj{t) (see 
fig lib)- The creep dynamics increases with the applied 
shear stress and decreases for larger values of the packing 
fraction. To quantify the creep dynamics, strain rates 
were computed for a time step St = Is and averaged 
over a time window of AT = 400s. On figEK, we dis- 
play 7 ^ — 7q ^ as a function of time for a fixed packing 
fraction at different applied shear stresses; 70 being the 
initial shear rate. This representation is a natural choice 
to probe a long-time logarithmic dynamics (consistent 
with 7 oc 1/t). Indeed, we observe a relation of the type: 
7 ^ — 7q ^ = Ct, corresponding to a long time logarithmic 

creep : 7 = 70 + ^ ln(l + Ct). On fig|3]3, we represent 

the values of C extracted from a linear fit as a function of 
G{(j))/a^ for different packing fraction values. We observe 



a monotonic increase, more prononced at larger packing 
fraction. 

Rheological model - To analyze quantitatively the data, 
we use a theoretical model introduced by Derec et al.jsl 
in the context of complex fluids rheology. We propose 
here to adapt it to the creep flows of granular packing. 
The model extends naturally the standard Maxwell visco- 
elastic rheology. It introduces an internal phenomeno- 
logical variable / called "fluidity" whose dimension is an 
inverse time. Physically, fluidity is a rate of stress relax- 
ation. To describe a complex dynamics displaying age- 
ing and rejuvenation, Derec et al. propose that fluidity 
should be time-dependent and its dynamics described by 
a simple " la Landau" phenomenological equation : 

dta = -fa + Gj (1) 
dtf = -af + n (2) 

The second equation introduces two dimensionless and 
positive parameters a and r. The first term corresponds 
a fluidity decrease i.e. an ageing process which renders 
the fluid more "viscous" with time. We call a the "age- 
ing parameter" . The second term, corresponds a fluidity 
increase (less "viscous") due to shear. We call r the "re- 
juvenation" parameter. The forms assumed by the age- 
ing and rejuvenation terms are the most simple one can 
get in an expansion consistent with a non-trivial long- 
time dynamics (see discussion in |6|). At constant shear 
stress cr, one obtains the relation fa = G7, yielding : 
dtf = -a(l — ( — )^)/^, where cfd = GJajv is the 

dynamical shear corresponding to steady shear rate and 
steady fluidity. Introducing an equivalent ageing param- 
eter : 

ae, = a(l-( — )2) (3) 

fo 

the solution of this equation is then : f(t) = — . 

1 + aeqft 

The shear rate variation is thus : 7 — 7o = ^cteq-t 
which leads to a long time logarithmic creep as observed 
experimentally. The experimental slopes G of fig. (jSb) 
can then be identifed using the relation : 

= eg (4) 

For all experiments performed at different stresses and 
packing fractions, the initial fluidity value fo = G{(j))jo/o' 
can be plotted as a function of (j). For all stresses, data 
cohapse onto a quasi-linear curve (see inset of figH]); the 
denser is the packing, the less is the initial fluidity. The 
hnear extrapolation of this curve to /o = yields a value 
(j)o = 0.635 ± 0.002, close to random close packing of 
monodisperse spheres. This can be interpreted as an 
arrest of the creep dynamics at a packing fraction cor- 
responding to the jamming limit for a random assembly 
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of frictionless spheres [l6|. Furthermore, the ageing dy- 
namics can be characterized by computing the equiva- 
lent ageing parameter a^^q according to relation (H]). If 
a and r are independent of shear, a^q should decrease 
quadratically and reach a zero value at a finite stress 
corresponding to the dynamical stress gd according to 
relation ([3]). Interestingly, the initial fluidity /o sets a 
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FIG. 4: Rescaled ageing parameter aeqfo as function of 
rescaled of shear stress dotted line: y — ^-(l — (i^)^) 
with t* = 250s. Inset : initial fluidity /o as a function of pack- 
ing fraction for different shear stresses, the color index reflects 
the ratio a /gm^ fit line: y — F{(j)Q — (/)) with F — 0.0086 and 
00 = 0.635 ± 0.002. 

time scale which can also be interpreted as an effective 
viscosity t^o = fo{(j)) diverging when approaching 

the jamming threshold from below : r]o oc {(j)o — 
Since, in principle, we are well below any thermalized 
regime where the viscosity concept could apply, this re- 
sult is quite remarkable. Moreover, we can use this orig- 
inal time scale to rescale the effective ageing parameter 
and obtain a collapse of all the data. On fig.(|4]), we plot 
Y = foaeq as a function of the non-dimensionalized stress 
: X = a/ar. The striking feature is that all data collapse 
onto a single curve for the entire range of stresses and 
packing fractions studied. The second important feature 
is that aeq displays a non monotonic behavior with a min- 
imum value at <j ~ (Tr (X = 1) corresponding to the on- 
set of the strain rate bursts identified in the stress-ramp 
experiments. This behavior is clearly the signature of in- 
ternal granular reorganizations leading to a slowing down 
of the creep dynamics instead of an increase as one might 
expect when shear is increased. For values above X ^ 2, 
the creep dynamics increases again before reaching the 
dynamical stress threshold {aeq = 0) at (Jd ^ 2.4(7^. For 
the smaller values of shear stress, i.e. below dr, the pre- 
dictions of Derec's model with constant coefficients can 
still be validated with a dynamical stress = 1.5(7^ (see 
dashed line on figH]). 

This experimental study shows that down to vanishing 
low applied shear and up to the yield stress value, inter- 
nal relaxation processes are present in a granular packing. 
The logarithmic creep hence observed, was analyzed us- 



ing a simple visco-elastic model which introduces a time 
dependent rate of relaxation (the fluidity). The dynamics 
is viewed as a competition beween intrinsic ageing and 
shear stress rejuvenation. The model allows a dynami- 
cal characterization of the initial packing fluidity which 
decreases linearly with packing fraction and vanishes at 
the random close packing limit. Under finite stress, we 
identified the onset of internal reorganizations, slowing 
down the creep process and setting the yield stress to 
higher values. This process could be related to the onset 
of shear induced anisotropy [17] or shear band formation. 
This internal dynamics is the sign of a peculiar fragility 
of this type of solid, possibly mediated by thermal acti- 
vation or by background mechanical noise. It may also 
be related to the intrinsic nature of the plastic response 
in amorphous solids [3], which are sometimes described 
as structurally fragile under finite shear [20]. 
We acknowledge the financial of the support CNRS- 
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